Novel structures of Co-Cu bimetallic clusters and their magnetic properties 
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The structural and magnetic properties of Coi8-mCu m 
(0 < m < 18) clusters are investigated with a genetic al- 
gorithm and a spd-band model Hamiltonian in the unre- 
stricted Hartree-Hock approximation respectively. In general, 
Cu atoms tend to occupy the surface, while Co atoms prefer 
to the interior of the clusters. Layered structures appear in 
some clusters with given stoichiometric compositions. The in- 
troduction of Cu atoms leads to large increase of the magnetic 
moment of Co-rich circumstance and nearly zero magnetism 
of the Cu-rich ambient. The interaction between Cu and Co 
atoms induces nonzero magnetic moment for Cu atoms. The 
total magnetic moments tend to decrease with the increase of 
Cu atoms. However, some particular large magnetic moment 
are found to be closely related to the structures. The envi- 
ronment of Cu and Co atoms have a dominant effect on the 
magnetism of the cluster. 
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Bimetallic clusters are an exciting research field due 
to their potential applicationsJjO the automobile indus- 
try and oil refined as catalystsau. Such nanoscale alloys 
may present a number of structures and phases that are 
different from those of corresponding pure metals. Previ- 
ously, there are intensive studies on homogeneous metal- 
lic clusters. But the reports on the bimetallic clusters are 
scarcest], especially for transition-metal bimetallic clus- 
ters because of the complexity in their electronic struc- 
ture. 

In this report, we exploit the structural and magnetic 
properties of bimetallic Co-Cu clusters. The main reason 
for choosing Co-Cu is that the physical properties of bulk 
Co and Cu are very different. We may get a clear picture 
of the various properties of the bimetallic clusters versus 
the different composition ratios. Moreover, Co and Cu 
alloys are non-miscible. The clusters may give a qualita- 
tive analysis from the mecoscopical points. 

Although the reliable results on clusters can be ob- 
tainedpOn. basis of quantum chemistry or density function 
theorytTcl, the well-known NP problem leads to expen- 
sive computational costs. Alternatively, empirical po- 
tential fitted from the bulk materials have been exten- 
sively employed .-tp study the structures and properties 
of clustersQU'EHrtj. In this letter, we obtain the lowest 
energy structures of Coi8- m Cu m (0 < m < 18) by a 
genetic algorithm (GA) with a Gupta-like many-body 
potentialtll. The parameters for inhomogeneous Cu-Co 
interaction are derived from the average pf. the Cu-Cu 
and Co-Co parameters. In the GA schemaljll3, a num- 



ber of random initial configurations are generated in the 
beginning. Then any two candidates in the population 
can be chosen as parents to generate a child cluster by 
mating operation. The obtained child cluster can be se- 
lected to replace its parent, if it has lower binding energy 
but its configuration is different from any one in the pop- 
ulation. 

We have checked the validity of current parameteriza- 
tion by ab initio calculation on the smallest clusters, i.e., 
homogeneous and inhomogeneous dimers and trimers. 
The ab initio calculation is performed by using DMqI 
package based on density functional theory (DFT)E3. 
During the DMol electronic structure calculations, the ef- 
fective core potential (ECP) and a double numerical basis 
including d-polarization function (DND) are chosen. The 
density function is treated Jijithin the generalized gra- 
dient approximation (GGA)Ej with exchange-eorrelation 
potential parameterized by Wang and Perdewc^l. A direct 
comparison of the ab initio and empirical results on the 
structural information such as equilibrium bond length 
and bond angle for those small clusters is given in Table. 
I. One can find that all the bond angles of either ho- 
mogeneous or inhomogeneous trimers are well described 
by empirical potential. Except C02 dimer, the difference 
of bond length between DFT and empirical calculation 
is less than 0.05A. We have further verified our empir- 
ical scheme by calculating the clusters C013 and Cui 3 . 
The average bond length and average binding energy pes 
atom of Co are 2j45A and 3.22 eV, which are 2.44 AEl 
and 3.66 ±0.36 eWu from first-principles calculations. For 
CU13, the bond length and the binding energy are 2.50A 
and 2.59 eV , in agreement with the results of TB-LMTO, 
2.52A and 2.46 eVEl. From the above comparisons, the 
overall agreement of our present model potential with ac- 
curate ab initio is rather reasonable. Therefore, we can 
use such Gupta-like potential in the global structural op- 
timization of 18-atom Cu-Co bimetallic clusters, in which 
ab initio calculations up to long time scale is computa- 
tional prohibitive. 

Table I. Bond length (d) and bond angle (8) are compared 
with the spin-polarized DFT-GGA method (in parenthe- 
sis) for small Co, Cu or Co-Cu clusters. 



Co 2 



Cu 2 



CuCo 



d (A) 2.203 (2.107) 2.235 (2.225) 2.219 (2.258) 

C03 Cu 3 C02CU CU2C0 

d (A) 2.230 (2.262) 2.337 (2.348) 2.308 (2.361) 2.324 (2.388) 

6 60.0 (60.0) 60.0 (60.0) 61.4 (62.1) 58.8 (56.5) 
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Fig.l gives the morphology structures corresponding 
to different stoichiometric composition of Cois-m(0 < 
m < 18). Great modifications are found in the bimetal- 
lic clusters. The most stable structures for the clusters 
with m = 0, 2, 15 — 18 are double icosahedron minus an 
atom in the layer, while the rest ones prefer to the bell- 
like structures. These imply that the mixing process has 
a great influence on the ground state structures. This 
may be original from the fact that 18-atom constitutes 
a double icosahedron with a defect. The existence of a 
defect may easily induce a structural transition. Hence, 
we can alter the composition ratios to attain some new 
structures. 




m=15 m=16 m=17 

FIG. 1. Optimized structures of Cois- 
concentrations. 
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Segregation effect is found in the ground state struc- 
tures of bimetallic clusters: Cu atoms tend to occupy 
the surface, while Co atoms prefer to the interior. The 
early occupied sites by Cu atoms have lower coordina- 
tion number (CN), then the higher CN. Take C05CU13 
as an example, eleven Cu atoms occupy the sites with 
CNs being 5 or 6, the rest 3 Cu atoms are in the layer 
with CNs=8. Moreover, the same kinds of atoms tend 
to assemble together in the same layer. The assembling 
Cu atoms tend to maximize the number of Cu-Cu and 
Co- Co bonds. This phenomenon may be due to the great 
difference of the surface energy and the cohesive energy 
between Co and Cu. To minimize the total energy, the 
atom with the smaller surface energy and cohesive en- 
ergy tends to occupy the surface, while the atom with a 
higher surface energy and cohesive energy favors to the 
interior. The average cohesive energy and surface energy 



of the bulk Cu, 3.544 eV and 1.934 Jm 2 , are smaller 
than those of the bulk Co, 4.386 eV and 2.709 Jnr 2 . 
Another possible reason is the atomic size effect. In our 
simulations, the first-nearest distance of Cu is 2. 556 A , 
larger than that of Co 2.507A. Thus, Co are more easily 
surrounded by Cu atoms. 

Another interesting finding is the appearance of lay- 
ered structures in the bimetallic clusters. For C013CU5, 
five Cu atoms occupy the lower CN sites while thirteen 
Co atoms constitute an icosahedron. Similar features are 
found in the clusters Co 8 Cui and Co7Cun,etc. These 
imply the existence of some ordering effects to maximize 
the number of Co-Cu bond. The ordering effect and the 
segregation effect coexist and compete with each other 
in the cluster, thus lead to the appearance of layered 
structure and the segregation of Cu atoms. The bizarre 
structure characters may bring some peculiar properties. 
Our previous studies have shown some peculiar thermal 
behavior in Cu-Co bimetallic clusterstfl. In the following, 
we investigate the electronic and magnetic properties by 
parameterized unrestricted Hartee-Fock approximation. 

The Hamiltonian, written in a local orbital basic set, 
has the expression: 



ji AT A 



(1) 



where c iacr { Cip a ) are the creation (annihilation) op- 
erators and n iaa is the number operator of an elec- 
tron.. The t"j is the hopping integral between different 
sites and orbitals. The orbital state a involved in cal- 
culation includes s,p x ,p Vl p Zl d X y 1 dy Zl d xz ,d x 2^ y 2,d :iz 2_ r 2. 
The single-site energy ti aa is given by 



e iaa = 4 + UAn(i) - \<rJKi) + £ An 0')Vy (2) 

Here e d refers to the orbital energy levels in the param- 
agnetic solutions of the bulk. An(j) denotes the charge 
change. The Coulomb interaction is described as 



Vu = 



U 



l + (C/i?y/e 2 ) 



(3) 



The orbital energy and the hopping integrals are taken 
to be the bulk values obtained from Andersen's linear 
muffin-tin orbital atomic sphere approximation(LMTO- 
ASA) paramagnetic bandstJ. The hopping are assumed 
to be spin independent and are averaged for the het- 
eronuclear. Exchange integralS|-ejther than J dd are ne- 



glected and JddiCo) = 0.99& 
U dd (Co) is obtained from Ref.f 
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The direct integral 
and U, 



/U dd relations 



are from the atomic tables. We take U, 



and U„ 



sd pd r 

come from Ref.t 
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s 1, Udd c3. For Cu, all the parameters 
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The magnetic moment can be determined by inte- 
grating the majority and minority local densities of 
state(LDOS) up to Fermi energy: 



Mi 



[p MT (e) - Pt ai (e)]de, 



(4) 



The LDOS is directly related to the diagonal elements 
of the local Green function by means of the recursion 
methocO: 
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Fig. 2 gives the total magnetic moment of Coi 8 - m Cu m 
bimetallic clusters as a function of the concentration of 
Cu atoms. The total magnetic moments decrease from 
33.61/if, to Ofib as the Cu concentrations vary from m = 2 
to m = 13. The curve can be divided to three sec- 
tions. The first section is the Co-rich circumstance with 
to = 1 — 5, where the introduction of a small amount of 
Cu atoms enhances the magnetism of the clusters. The 
contribution to the magnetism mainly comes from the 
Co atoms far away from Cu atoms. It may be due to 
the large charge transfer between Cu and Co atoms. The 
second section is the comparative composition ratio with 
to = 6 — 12, in which the magnetic moments fluctuate 
with the cluster size. For the clusters with concentra- 
tions to = 6 — 8, the case is very similar to the first. 
But in the case of to = 9 — 11, the interior pentagonal 
bipyramid significantly contributes to the magnetism of 
the cluster. Thus, relatively large magnetic moments are 
found for these clusters. The third section is the Cu- 
rich clusters with to = 13 — 17. Cu atoms seem to have 
a "screen" effect on the magnetic moments and lead to 
zero magnetic moment in these clusters. 
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FIG. 2. The total magnetic moments of Coi8- m Cu m 
bimetallic clusters as a function of the concentrations of Cu 
atoms. 

The magnetism of the cluster is found closely related 
to the environments of Co and Cu atoms and the cluster 
geometrical characters. For example, the total magnetic 
moment of C016CU2 is much larger than those of C017CU1 
and C015CU3. The main reason is due to the large charge 



transfer between Co-Cu. Obviously, the number of Co- 
Cu bond in C016CU2 is much more than that of C017CU1 
and C015CU3. The more the number of Co-Cu bond in 
the cluster, the larger the charge transfer, which induces 
a large magnetic moment. Similarly, the large magnetic 
moments for C011CU7 and CosCuio are obtained, com- 
pared with CoioCus and C07CU11. Contrary, the to- 
tal magnetic moments for C013CU5 is found much larger 
than that for C012CU6, even the number of Co-Cu bond 
in the former is more than the latter one. As discussed 
above, Coi 3 Cu 5 is a good layered structure with 5 Cu 
atoms in the same layer and 13 Co atoms constituting 
an icosahedron. Although the high symmetry compresses 
the magnetism of the clusters, the 5 Cu atoms has an 
less effect to the icosahedron. It is well-known that the 
smaller the cluster size, the larger the magnetic moments. 
Therefore, a large magnetic moment is obtained for the 
cluster C013CU5. 
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FIG. 3. Average magnetic moments of Co and Cu atoms 
as a function of the concentrations of Cu atoms. 



It is also worthy to note the nonzero magnetic moment 
of Cu atoms in some bimetallic clusters. Fig. 3 shows 
the average moments of Co and Cu as a function of the 
concentration of Cu. Although the average magnetic mo- 
ment of Cu atoms is zero or nearly zero in most cases, par- 
ticularly high magnetic moments are found in the clusters 
with to = 1 — 3 and less pronounced peaks are found at 
to = 6, 12. For the case of Co atoms, the hybridization 
with the Cu atoms leads to an oscillatory behavior for 
the average magnetic moment. These may be due to the 
different charge transfer. For Co atoms, the charges are 
transferred from sp orbitals to d orbital, contrary to the 
case of Cu atoms that the charges are transferred from d 
orbital to sp orbitals. Further, the charge transfer takes 
place from Co atoms to Cu atoms for the clusters with 
Cu concentration to < 13, while the case is reverse for 
to > 13. For the case of to = 3, 6, 12, the charge transfer 
from Co atoms to Cu atoms is found very large, which 
induces the large magnetism of Cu atom. 

To explore the origin of peculiar magnetic properties, 
we show the total density of states(DOS) and sp,d DOS of 
the pure and represented bimetallic clusters in Fig. 4. The 
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cluster Fermi level is presented as a dashed vertical line 
and shifted to zero. In general, the DOS near to Fermi 
level play a primary role in determining the magnetism of 
the clusters. Obviously, the contribution of d electrons is 
found dominantly, while the sp electrons contribution is 
low. The contribution of d electrons in C016CU2 is larger 
than that in Cois, which leads to a large increase of the 
magnetic moments in C016CU2. Similarly, the contribu- 
tion of d electrons in C016CU2 is also larger than that 
in C013CU5, thus their corresponding magnetism are dif- 
ferent from each other. For Co 10 Cu 8 , the contribution 
of d electrons is relatively less near to the Fermi level 
compared with other cases, which leads to a particular 
small magnetic moments. Moreover, the hybridization 
between sp and d orbitals among Co-Cu, Co-Co and Cu- 
Cu atoms in CoioCus is also smallest among these four 
cases, while it is largest in C016CU2. This also enhances 
the magnetism of C016CU2. 




Energy (eV) Energy (eV) 

FIG. 4. Total, sp and d orbitals DOS for some compositions: 
(a) Cois; (b) Coi6Cu2;(c) C013CU5, (d)CoioCus; The vertical 
dashed lines indicate the DOS integrated Fermi level. 



In summary, the geometrical and magnetic properties 
of bimetallic clusters Coi 8 _ m Cu m have been studied by 
a genetic algorithm and a spin polarized tight-binding 
Hamiltonian. The main conclusions can be made as fol- 
lows: (1) Great modifications are found for the 18-atom 
bimetallic clusters due to different composition ratios. 
This suggests that we can alter the composition ratios 
to attain new structures. (2) Ordering effect and seg- 
regation effect influence the atomic configurations of the 
bimetallic clusters simultaneously, which leads to the seg- 
regation of the Cu atoms to the surface and some layered 
structures. These may explain why there are no corre- 
sponding ordered compounds of Co-Cu bulk in the low 
temperature. (3) The introduction of Cu atoms causes a 
dramatic increase of magnetism in Co-rich circumstance 
and nonzero moments for Cu atoms. Particular large 
moments can be associated with the environments of Co 
and Cu atoms and the geometrical characters. (4) Cu 
atoms have an "screen" effect on the cluster magnetism 
in Cu-rich ambient and enhance the magnetics in Co-rich 



environment. 
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